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ABSTRACT 

The most metal-poor DLA known to date, at z a b s = 2.61843 in the spectrum of 
the QSO Q0913+072, with an oxygen abundance only ~ 1/250 of the solar value, 
shows six well resolved Di Lyman series transitions in high quality echelle spectra 
recently obtained with the ESO VLT. We deduce a value of the deuterium abun- 
dance log(D/H) = —4.56 ± 0.04 which is in good agreement with four out of the six 
most reliable previous determinations of this ratio in QSO absorbers. We find plausi- 
ble reasons why in the other two cases the la errors may have been underestimated 
by about a factor of two. The addition of this latest data point does not change 
significantly the mean value of the primordial abundance of deuterium, suggesting 
that we are now converging to a reliable measure of this quantity. We conclude that 
(log(D/H) p ) = -4.55 ±0.03 and f2 b ,o^ 2 (BBN) = 0.0213 ± 0.0010 (68% confidence 
limits). Including the latter as a prior in the analysis of the five year data of WMAP 
leads to a revised best-fitting value of the power-law index of primordial fluctuations 
n s = 0.956 ±0.013 (la) and n s < 0.990 with 99% confidence. Considering together the 
constraints provided by WMAP 5, (D/H) p , baryon oscillations in the galaxy distri- 
bution, and distances to Type la supernovae, we arrive at the current best estimates 
n hfi h 2 = 0.0224 ± 0.0005 and n s = 0.959 ± 0.013. 
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1 INTRODUCTION 

The exquisite precision with which the temperature 
anisotropies of the cosmic microwave background (CMB) 
have been mapped on the sky has allowed the determina- 
tion of cosmological parameters to better than 10% (Dunk- 
ley et al. 2008). It is still important, however, to measure 
these parameters with alternative methods, partly as a con- 
sistency check on the standard cosmological model, but also 
because the CMB fluctuations generally constrain combina- 
tions of more than one parameter (e.g. Bridle et al. 2003). 
The contribution of baryons to the present-day critical den- 
sity, where as usual h — J/o/100km s _1 Mpc -1 , is 
one such example: the value of Qh,oh 2 which best fits the 
power spectrum of CMB fluctuations is tied to other pa- 
rameters, such as the spectral index of primordial density 
perturbations, n s , and the optical depth to reionisation, r 
(see, for example, Pettini 2006). 

The quantity Q^^k 2 can also be deduced from the pri- 
mordial abundances of the light elements whose nucleosyn- 



thesis is the result of physical processes that are entirely 
different from the acoustic oscillations of the photon-baryon 
fluid imprinted on the CMB, and that took place at much 
earlier times — only a few hundred seconds, rather than a 
few hundred thousand years, after the Big Bang. Among 
the light elements produced by primordial nucleosynthesis, 
deuterium is the one whose abundance by number relative to 
hydrogen, (D/H) p , depends most sensitively on fl bt0 h 2 , and 
is thus the 'baryometer' of choice (Steigman 2007; Molaro 
2008). 

Reliable measures of (D/H) p are difficult to obtain, 
however. The astrophysical environments which seem most 
appropriate are the hydrogen-rich clouds absorbing the light 
of background QSOs at high redshifts, but rare combina- 
tions of: (a) neutral hydrogen column density in the range 
17 < log[iV(Hl)/cm~ 2 ] < 21; (b) low metallicity [M/H] cor- 
responding to negligible astration of D; and (c) most im- 
portantly, low internal velocity dispersion of the absorbing 
atoms allowing the isotope shift of only 81.6 km s" 1 to be 
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adequately resolved, are required for this observational test 
to succeed. Thus, while the potential of this method was 
appreciated more than thirty years ago (Adams 1976) and 
first realised with the advent of 8-10 m class telescopes in 
the 1990s (Tytler et al. 1995), the number of trustworthy 
measurements of (D/H) p is still only five or six (O'Meara et 
al. 2006; Steigman 2007). 

There are a strong incentives to increase these mea- 
gre statistics. On the one hand, the mean ((D/H) p ) = 
(2.82 ±0.26) x 10 -5 implies « b , /!. 2 (BBN) = 0.0213 ± 0.0013 
(O'Meara et al. 2006) which agrees, within the errors, with 
the value fi b ,oh 2 (CMB) = 0.02273 ± 0.00062 obtained by 
Dunkley et al. (2008) from the analysis of five years of ob- 
servations with the Wilkinson Microwave Anisotropy Probe 
(WMAP 5). On the other hand, the standard deviation 
among the six measures of (D/H) p in high redshift QSO 
absorbers exceeds the dispersion expected from the individ- 
ual errors (leading to a high value of \ 2 of t ne srx measure- 
ments about the weighted mean). This is probably due to 
the considerable difficulties in accounting for the full mea- 
surement errors (random and systematic) of astronomical 
observations, but is nevertheless a source of concern, par- 
ticularly when viewed in conjunction with the as yet poorly 
understood dispersion of D /H values in the Milky Way (Lin- 
sky et al. 2006). 

Additional measurements of the D /H ratio at high red- 
shift should, in principle at least, lead to a more precise 
estimate of (D/H) p , as well as to a better assessment of the 
reasons for the scatter of the existing values. Such improve- 
ments are not only of interest in a cosmological context, but 
are also relevant to Galactic chemical evolution models in 
which the degree of astration of D is an important diagnostic 
(e.g. Romano et al. 2006; Steigman, Romano & Tosi 2007). 
In this paper, we report observations of deuterium absorp- 
tion in a seventh high redshift absorber, a very metal-poor 
damped Lya system (DLA) at z abs ~ 2.61843, and consider 
their implications for the determination of (D/H) p , Qb,oh 2 , 
and other cosmological parameters^ 



2 OBSERVATIONS 

The existence of a metal-poor DLA at z a hs = 2.61843 in the 
line of sight to the bright (V = 17.1), z cm = 2.785, QSO 
Q0913+072 has been known for some time (Pettini et al. 
1997; Ledoux et al. 1998; Erni et al. 2006). More recently, 
it has been realised that DLAs of low metallicity are also 
likely to have simple kinematics, possibly reflecting an un- 
derlying mass-metallicity relation of the host dark-matter 
halos (Ledoux et al. 2006; Murphy et al. 2007; Prochaska et 
al. 2008; Pontzen et al. 2008). The most metal-poor DLAs 
are thus likely to be among the most promising candidates 
for follow-up high-resolution spectroscopy aimed at resolv- 
ing the isotope shift in high order lines of the Lyman series. 

To this end, we targeted Q0913+072 with a concerted 
series of observations in 2007, using UVES (Dekker et al. 
2000) on the ESO VLT2. The data were acquired in service 

1 Damped Lya systems are a class of QSO absorbers de- 
fined by their high column densities of neutral hydrogen, 
log[iV(Hl)/cm- 2 ] ^ 20.3; see Wolfe, Gawiser, & Prochaska 
(2005), and Pontzen et al. (2008) for recent reviews of DLAs. 



mode; including a few spectra retrieved from the UVES data 
archive, the total exposure time devoted to this QSO was 
77 550 s. Details of the data reduction process can be found 
in Pettini et al. (2008); briefly, the final reduced and com- 
bined spectrum covers the wavelength region from 3310 A 
to 9280 A with a resolution of 6.7 km s -1 FWHM, sampled 
with ~ 3 pixels. The signal-to-noise ratio in the continuum 
varies from S /N ~ 35 per pixel at the redshifted wavelength 
of the damped Lya line, to S/N ~ 15 just longwards of the 
Lyman limit of the z a b s = 2.61843 DLA. 

Pettini et al. (2008) presented an analysis of the chemi- 
cal composition and kinematics of this and other metal-poor 
DLAs, from consideration of metal absorption lines due to 
H, C, N, O, Al, Si, and Fe; we refer the interested reader 
to that paper for details. Of relevance to the present work 
are the following main results. With relative abundances 
[C/H] = -2.75 and [Fe/H] = -2.80 the z abs = 2.61843 sys- 
tem in Q0913+072 is the most metal-poor DLA known0 Its 
oxygen abundance, [O/H] = —2.37, is also among the lowest 
recorded in damped systems. The metal lines arising in the 
neutral gas have a very simple kinematic structure, consist- 
ing of two components at z a b s = 2.61828 and z a b s = 2.61843 
(Av = 12.4km s _1 ), with low internal velocity dispersions, 
b — 3.7 km s _1 and 5.4 km s" 1 respectively (where b = ~j2a 
and a is the one dimensional velocity dispersion of the ab- 
sorbing atoms along the line of sight). The column density 
of neutral gas is distributed between the two components in 
approximately 3:7 proportion (see Figure 1 and Table 2 of 
Pettini et al. 2008). 

While these characteristics bode well for the possibility 
of resolving the D component of the Lyman series lines, they 
do not guarantee it. It is often the case that partially ionised 
gas at nearby velocities, with optical depths too low to be 
recognised in the metal lines, can blend with the H and D 
absorption due to the DLA itself, and prevent a reliable 
measurement of 7V(Dl). Partially ionised gas does in fact 
appear to contribute to the C II absorption lines in this DLA, 
but seems to be confined to positive, rather than negative, 
velocities relative to z abs = 2.61843. As can be seen from 
Fig. [1] this is indeed, somewhat fortuitously, the case. 

3 HYDROGEN AND DEUTERIUM 

ABSORPTION IN THE Zabs = 2.61843 DLA 

3.1 Deuterium 

In Fig. [T]we have reproduced on a velocity scale relative to 
2 abs = 2.61843 the normalised profiles of absorption lines 
in the Lyman series of the DLA, from Lya to LylO. Lyll is 
also visible in the LylO panel, while higher order lines are all 
blended with one another, constituting an 'effective' Lyman 
limit from Ao = 917 A. It is clear from the high order transi- 
tions that there is Hi absorption at positive velocities rela- 
tive the z a bs = 2.61843 DLA, extending to v ~ +200 km s^ 1 . 
On the blue side, however, there appears to be no additional 
components over those seen in the metal lines, and the D I 
absorption is clearly resolved in Ly<5, and Ly7 through to 
Lyll, six transitions in total. (In Lya, Ly/3 and Ly7, which 



2 With the usual definition, [X/H] = log(X/H) DLA - log(X/H) . 
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Figure 1. Observed profiles (black histograms) and fitted Voigt profiles (continuous red lines) of absorption lines in the Lyman series 
of the 2 a t, s = 2.61843 DLA in Q0913+072. The j/-axes of the plots show relative intensity. The two vertical tick marks in each panel 
indicate the expected locations of the main absorption component of the DLA in, respectively, Hi (at v = 0km s — 1 in the plots) and 
Dl (at v = —81.6km s — 1 ). A second absorption component, centred at v = —12.4km s — , is resolved in the metal lines associated with 
this DLA, but the two components are blended in the intrinsically broader H I and D I absorption. Additional H I absorption is seen 
at positive velocities relative to the DLA, but its column density is ~ 100 times lower than that of the DLA, and therefore does not 
contribute to the observed Dl absorption lines (see text). 
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Table 1. Absorption components fitted to the z a b s = 2.61843 
DLA in Q0913+072. 



ion 




b 

(km s" 1 ) 


Fraction 


Hi a 


2.61828 b 


11. 9 C 


0.33 d 


Hi a 


2.61843 b 


11.1° 


0.67 d 


Di 


2.61828 b 


8.6 


0.33 


Dl 


2.61843 b 


8.7 


0.67 


Oi 


2.61828 


3.5 


0.31 


Oi 


2.61843 


5.2 


0.69 



a Additional components at positive velocities relative 
to z a bs = 2.61843 were fitted to the Hi lines. 

b Fixed to be the same as Oi. 

c Poorly determined because all Hi lines are 
saturated and blended. 

d Fixed to be the same as D I. 

have the highest transition probabilities, Hi and Dl are in- 
trinsically blended, whereas in Lye and Ly6 presumably un- 
related absorption in the Lyman forest is blended with the 
absorption due to the DLA). 

The availability of many D I transitions with a range of 
optical depths, from the saturated Ly<5 to the optically thin 
LylO and Lyll, allows for a precise determination of the 
column density N(D i). To this end, we fitted the profiles of 
the Lyman lines where the D I absorption is resolved with 
theoretical Voigt profiles generated by the VPFIT (version 
8.02) software package0 VPFIT uses \ 2 minimisation to 
deduce the values of redshift z, logarithmic column density 
log[7V/ (cm -2 )], and Doppler parameter b (km s _1 ) that best 
reproduce the observed absorption line profiles, taking into 
account the instrumental broadening function in its % 2 min- 
imisation and error evaluation. We used the compilation of 
laboratory wavelengths Ai a b and /-values by Morton (2003). 
In Table [T] we have collected relevant details of the absorp- 
tion model fitted to the observed line profiles in Fig. [T] the 
corresponding theoretical profiles generated by VPFIT are 
shown with continuous red lines in Fig. [T] 

As mentioned in Section 2, the metal lines detected in 
this DLA consist of two components (Pettini et al. 2008), 
with the parameters listed in Table Q] for Oi. The two com- 
ponents are not clearly distinguished in D I because: (i) ab- 
sorption by the lighter D is intrinsically broader (see discus- 
sion below), and (ii) the S/N ratio of the spectrum is lower 
at the shorter wavelengths of the high order Lyman lines (at 
z a b s = 2.61843 the wavelength region between Ly<5 and the 
Lyman limit is redshifted to A b s = 3438-3318 A, approach- 
ing the cut-off by atmospheric ozone). Thus, in the first step 
of the fitting process, we fixed the redshifts of the two com- 
ponents to be the same as those determined from the metal 
lines while keeping the values of b and A^(D i) as free param- 
eters to be optimised by VPFIT. The best fitting values of b 
and N(Dl) thus reached are consistent with those deduced 
for O I and other metal lines. The relative proportions of 
atoms between the two components are similar in all species 
and the b- values are broader in D than in O, as expected. Re- 
call that the velocity dispersion of an absorption line is the 

3 VPFIT is available from [http: / /www.ast.c am.ac.uk/~| 
rfc/vpfit.html 



quadratic combination of macroscopic (turbulence) and mi- 
croscopic (temperature related) terms, i.e. b 2 ot = b^ UTb + 6^; 
while the former is presumably the same for all elements, 
the latter has an inverse dependence on the square root of 
the mass, since b\ = IkTjra (e.g. Stromgren 1948). The 
higher b- values of the two absorption components in D than 
O imply temperatures of a few thousand degrees. 

The model parameters listed in Table [T] correspond to 
a best fitting value for the column density of D I (sum of the 
two components) of 

log[AT(Di)/cm- 2 ] = 15.78 ±0.02, 

but we stress than this value does not depend on the de- 
tails of the 'cloud' model adopted. To test its robustness, 
we rerun VPFIT without any assumptions about the dis- 
tribution of velocities of the absorbers. Using as a start- 
ing point a single absorbing component with unspecified 
redshift and 6-value, VPFIT returned a best-fitting value, 
log[A r (D i)/cm~ 2 ] = 15.79±0.02, that is only 0.01 dex higher 
than that obtained with the parameters in Table[T] (although 
this second model results in a higher value of x 2 between 
computed and observed profiles). Any other adjustment to 
the details of the profile fits, such as varying the relative 
proportions of D I between the two components and their 
6-values, or changing the parameters of nearby absorption 
in the Lyman forest, resulted in even smaller differences in 
the total value of log JV(Di) returned by VPFIT. 

As a final comment here, we point out that the VPFIT- 
generated line profiles shown in Figure [1] include a number 
of other components to the red of z a b s = 2.61843, required 
to reproduce the saturated profiles of the Hi Lyman lines 
which extend to v ~ 200 km s" 1 . The parameters of these 
extra absorptions are poorly constrained, but do not influ- 
ence the determination of log N(D i) because they involve 
column densities of less than ~ 1/100 that of the DLA and 
thus do not contribute to the observed D I absorption. 

3.1.1 Consistency check with the apparent optical depth 
method 

Since the D I absorption lines we observe are resolved, we 
can deduce the column density directly from the measured 
residual intensity in each wavelength, or velocity, interval 
across the absorption lines (e.g. Hobbs 1974). The apparent 
optical depth at velocity v, T a (v), is related to the observed 
intensity in the line, I bs(v), by 

Ta (v) = -In [I obs (v)/I (v)l (1) 

where Iq(v) is the intensity in the continuum. With the as- 
sumption of negligible smearing of the intrinsic line profile 
by the instrumental broadening function, we have: 

r(v) » r a (v). (2) 

The optical depth t(v) is in turn related to the column 
density of D I atoms in each velocity bin, N(v) in units of 
cm -2 (km s -1 ) -1 , by the expression 

where the symbols /, A, c, e and m e have their usual mean- 
ings. 
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Table 2. Results of apparent optical depth analysis of D I lines 1 



Transition 



K K 

(A) 



7V tot (Di) 
(10 15 cm- 2 ) 



Dl Ly7 
Dl Ly8 
Dl Ly9 
Dl LylO 


925.9737 
922.899 
920.712 
919.102 


0.003184 
0.002216 
0.001605 
0.001201 


5 64+ - 43 

6 S5+ 40 

D.OO_q 2g 

6.08l° ;S 

5 61 + 041 
°' Di -0.37 


Weighted Mean 






5 93+ 21 



a Morton (2003). 

As emphasized by Savage & Sembach (1991), the attrac- 
tion of the apparent optical depth method lies in the fact 
that no assumption has to be made concerning the veloc- 
ity distribution of the absorbers. Furthermore, this method 
provides a consistency check when two or more transitions 
arising from the same atomic energy level, but with differ- 
ent values of the product /A, are analyzed, as is the case 
here. The run of N(v) with v should be the same, within 
the errors in I b B (v), for all such lines. 

For Dl in Q0913+072, we can apply to apparent optical 
depth analysis to four unsaturated transitions, from Ly 7 to 
Ly 10; the other lines in the D I Lyman series being either 
partly blended or saturated (see Figure [TJ. For each transi- 
tion, we deduced a value of N(D i) by summing eq. (|3} over 
the n velocity bins which make up the absorption profile: 



A TO T = ^Ni(v) 



(4) 



From the known error in the value of Idbs(v) in each velocity 
bin, 5I h s (y), we calculated the error SNtot- 



SN TOT =^T5Ni(v) 



(■>) 



which is asymmetric about Atot because of the non-linear 
nature of eq. {TJ. 

Figure [2] shows the run of Ni(v) across the four Dl ab- 
sorption lines, while values of Atot and SNtot are listed 
in Table [2] As can be seen from Figure [2] and Table [2j 
the four D I absorption lines are in good mutual agree- 
ment, and the standard deviation between the four inde- 
pendent measures of iV(Dl) is la = 0.35 x 10 15 cm -2 , or 
~ 6% of the mean. The weighted mean of the four values 
of N(D i) returned by the apparent optical depth analy- 
sis is (N(D I)) = (5.93+^ig) x 10 15 cm- 2 , where the un- 
certainties quoted are the errors on the weighted mean. 
On a logarithmic scale log[AT(D l)/cm" 2 ] = 15.775 ± 0.014, 
in very good agreement (as expected) with the value of 
log[AT(Dl)/cm- 2 ] = 15.78 ±0.02 returned by VPFIT, which 
we retain in the subsequent analysis. 



3.2 Hydrogen 

Fig. [3] shows the region encompassing the damped Lya line. 
In principle, the column density of neutral hydrogen should 
be tightly constrained from the shape of the damping wings 
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Figure 2. Apparent optical depth analysis: column density as 
a function of velocity for the four unsaturated and unblended D I 
lines in the z abs = 2.61843 DLA in Q0913+072. The transitions 
shown in the legend are in increasing order of /A from the top. 
Values of TVtot(Di) are collected in Table [2] 



which extend over many hundreds of pixels (and are thus ex- 
tremely well sampled by the data). In practice, the limiting 
factors are the uncertainty in the determination of the con- 
tinuum level (to which the absorption is normalised) and the 
overlapping absorption from the narrower lines in the Lya 
forest. After numerous trials varying the continuum level 
and the weights given to different spectral intervals deemed 
to be free of overlapping absorption (see Kirkman et al. 2003 
for a more extensive discussion of the problem) we converged 
on the best fitting value: 



Iog[AT(Hi)/cm" 



20.34 ± 0.04; 



the corresponding theoretical damped profile is overplotted 
on the data in Fig. [3] 

It is difficult to estimate reliably the systematic error 
which may be affecting this determination. While consistent 
with the above estimate of log iV(Hl) from Lya, the higher 
order Lyman lines do not help to reduce the error further 
because the are all saturated and their equivalent widths 
include uncertain contributions from the lower column den- 
sity components to the red of the DLA. In future, it may 
be possible to improve on the accuracy with which values 
AT (Hi) can be deduced from the analysis of damped profiles 
buried within the Lya forest by developing more sophisti- 
cated statistical methods to deal with the overlapping ab- 
sorption, perhaps analogous to those used to subtract fore- 
ground Galactic emission from maps of the CMB. For the 
moment, we can perhaps obtain an indication of the magni- 
tude of the systematic uncertainty in JV(H i) by considering 
three previous estimates of this quantity, reported by differ- 
ent authors who used different spectrographs and telescopes, 
as follows: log [AT (H i)/cm~ 2 ] = 20.36 ± 0.08 (Pettini et al. 
1997); 20.2 ± 0.1 (Ledoux et al. 1998); 20.36 ± 0.05 (Erni 
et al. 2006). Of these, only the last one refers to a (small) 
subset of the data used here. From these results it appears 
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Figure 3. Top: Portion of the normalised spectrum of Q0913+072 encompassing the damped Lyct line at z a j, s = 2.61843 (black 
histogram). The continuous red line shows the theoretical profile for a neutral hydrogen column density log[A r (H l)/cm -2 ] = 20.34. 
Bottom: Expanded portion of the upper plot on a scale chosen to illustrate the sensitivity of the damping wings to the column density 
of neutral hydrogen. The short-dash green lines correspond to the ±0.04 limits to the best fitting value log[7V(H l)/cm~ 2 ] = 20.34. In 
both panels the y-axis is residual intensity. 



that our present estimate, log [AT (H l)/cm~ 2 ] = 20.34 ±0.04, 
is unlikely to be in error by more than the stated margin. 

From the values of AT(Di) and AT(Hi) deduced in sec- 
tions 3.1 and 3.2, we arrive at the determination of the deu- 
terium abundance in the z a b s = 2.61843 DLA in line to 
Q0913±072 of: 

log(D/H) = -4.56 ±0.04 
(where the errors have been combined in quadrature). 



4 THE PRIMORDIAL ABUNDANCE OF 
DEUTERIUM 

In Table [3] we have collected relevant measurements for all 
the high redshift QSO absorption systems where the isotope 
shift has been resolved in absorption lines of the Lyman se- 
ries. This prime sample of what are generally considered to 
be the most reliable measures of D/H at high z now con- 
sists of seven independent determinations. Other reports in 
the literature (e.g. Levshakov et al. 2002; Chrighton et al. 
2004), while still interesting, refer to spectra of less straight- 
forward interpretation because not all the D I components 
are resolved, so that the values of D/H deduced are more 



dependent on the precise description of the kinematics of 
the gas than is the case for the seven absorbers listed in Ta- 
ble [3] (see also the discussions of this point by Kirkman et 
al. 2003 and O'Meara et al. 2006). 

All the measurements of D/H in Table [3] should be rep- 
resentative of the primordial abundance (D/H) p , because in 
all seven cases the gas has undergone little chemical enrich- 
ment, as evidenced by the low abundance of O (and other 
heavy elements), less than 1/10 of solar (see column 5 of 
Table For comparison, the total degree of astration of D 
over the lifetime of the Milky Way amounts to only ,$ 20%, 
if one accepts the possibility that in the local interstellar 
medium some of the deuterium may be depleted onto dust 
grains (Linsky et al. 2006). Theoretically, galactic chemical 
evolution models (e.g. Prantzos & Ishimaru 2001; Romano 
et al. 2006) show negligible reduction in D/H from the pri- 
mordial value when the gas metallicities are as low as those 
of the seven QSO absorbers in Table[3] while observationally 
no trend is observed between D/H and O/H (Pettini 2006). 
Dust depletion of D is not expected to be an issue here, given 
the very low depletions of even the most refractory elements 
at these low levels of chemical enrichment (e.g. Akerman et 
al. 2005). 

We arrive at an estimate of log(D/H) p by averaging 
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Table 3. Prime Sample of D/H measurements in QSO Absorption Line Systems 
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QSO z cm z abs logJV(Hi) [0/H] b log (D/H) Rcf. a 

(cm -2 ) 



HS 0105+1619 


2.640 


2.53600 


19.42 ±0.01 


-1.70 


-4.60 ±0.04 




1 


Q0913+072 


2.785 


2.61843 


20.34 ± 0.04 


-2.37 


-4.56 ± 0.04 


2. 


3 


Q1009+299 


2.640 


2.50357 


17.39 ±0.06 


< -0.67 c 


-4.40 ±0.07 




4 


Q1243+307 


2.558 


2.52566 


19.73 ±0.04 


-2.76 


-4.62 ± 0.05 




5 


SDSS J155810.16-003120.0 


2.823 


2.70262 


20.67 ±0.05 


-1.47 


-4.48 ± 0.06 




6 


Q1937-101 


3.787 


3.57220 


17.86 ±0.02 


< -0.9 


-4.48 ± 0.04 




7 


Q2206-199 


2.559 


2.07624 


20.43 ±0.04 


-2.04 


-4.78 ± 0.09 
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References— 1: O'Meara et al. (2001); 2: Pettini et al. (2008); 3: This work; 4: Buries & Tytler 
(1998b); 5: Kirkman et al. (2003); 6: O'Meara et al. (2006); 7: Buries & Tytler (1998a); 8: Pettini 
& Bowen (2001). 

b Relative to the solar value log(O/H) + 12 = 8.66 (Asplund, Grevesse & Sauval 2005). 
c This is a very conservative upper limit on the metallicity. Buries & Tytler (1998b) estimate 
[Si/H] ~ —2.5 and [C/H] ~ —2.9 from photoionisation modelling. 



the individual measure of log (D/H) in Table[3]to obtain the 
weighted mean (log(D/H) p ) = —4.55 ± 0.02. However, the 
scatter of the points about this value is rather high for the 
quoted error bars assuming a Gaussian error model, giving 
X 2 = 19 which formally corresponds to a high probability 
that a further independent experiment would obtain a better 
fit to to the data, P(x 2 > 19) < 0.01. This suggests that the 
errors on the individual measures of D/H may have been 
underestimated - a well known and much discussed (e.g. 
Steigman 2007) 'problem' of the determination of (D/H) p 
from QSO spectra. An alternative, and more conservative, 
method of estimating the error on the weighted mean would 
attempt to account for the actual scatter in the data. One 
possibility is to consider the mean and standard deviation 
of a large number of weighted means generated by random 
sampling (with substitution) of the seven observed values 
(the bootstrap method; see Efron & Tibshirani 1993). In 
this way we obtain: 

(log (D/H) p ) = -4.55 ±0.03 

A more Bayesian approach could try to account for uncer- 
tainty in the error bars using the observed scatter in the 
data. For example, fitting a Gaussian model with the vari- 
ance on each point increased by a constant to a 2 + 8a 2 shows 
that the data prefer values of 8a > 0. Using the maximum 
likelihood value 5a ~ 0.07 gives a posterior constraint on 
(log(D/H) p ) consistent with the bootstrap estimate above, 
as does marginalizing over 8a with prior oc l/(0.02 2 +8a 2 ) 
(though the distribution is somewhat non- Gaussian). Using 
a model that multiplicatively increases the noise gives simi- 
lar results. 

The scatter in the reported determinations of D/H in 
QSO absorbers is illustrated in Figure [4] Out of the seven 
measurements, there are two which lie outside the confidence 
intervals of the mean by more than la: Q1009+299 and 
Q2206— 199. If we consider these two cases more closely we 
can find, in retrospect, plausible reasons why their la errors 
may have been underestimated. The partial Lyman limit 
system in line to Q1009+299 has the lowest column density 
among the sample, log[N (Hi) /cm~ 2 ] = 17.39, and is unique 
among the ones considered here in showing D I absorption in 
only one line, Lya, the column density of D I being too low 
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Figure 4. Measures of deuterium abundance in high rcdshift 
QSO absorbers. Only cases were the deuterium absorption is 
clearly resolved from nearby spectral features are shown here (see 
text). Blue circles denote systems observed from the ground with 
8-10 m telescopes and echelle spectrographs, while the red trian- 
gle refers to lower resolution observations made with the Hubble 
Space Telescope. Absorption systems discussed in the text are 
labelled with the name of the background QSO. The horizontal 
lines are drawn at the weighted mean value of log (D /H) and its 
error, as determined with the bootstrap method. 



to produce discernible absorption in higher order lines. The 
2 abs = 2.07624 DLA towards Q2206-199 is the lowest red- 
shift absorber in the sample, requiring space-borne observa- 
tions to record the high order lines where D I absorption can 
be resolved from H I. The HST-ST1S spectrum of this object 
published by Pettini & Bowen (2001) is of lower S/N ratio 
and resolution than the other six cases in Table[3]which were 
all obtained with ground-based 8-10 m telescopes and echelle 
spectrographs. If we arbitrarily double the la estimates of 
the errors in the determinations of D/H in Q1009±299 and 
Q2206— 199 reported in the original works, we find the same 
weighted mean (log(D/H) p ) = 4.55 ± 0.02 as before, but a 
much reduced \ 2 = H {P{x 2 > H) - 0.1). 
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In conclusion, the value of D /H we deduce here for the 
Zabs = 2.61843 DLA in the spectrum of Q0913+072 is in 
good agreement with four out of the previous six determi- 
nations generally considered to be the most reliable. Its in- 
clusion in the 'prime' sample, helps identify two outliers and 
plausible reasons for much of the scatter among the sample. 
The weighted mean (log(D/H) p ) is unchanged compared 
to the most recent previous estimate of this quantity by 
O'Meara et al. (2006), leading us to conclude that the true 
value of the primordial abundance of deuterium lies in the 
range: 

(log (D/H) p ) = -4.55 ± 0.03 (6) 
at the 68% confidence level. 

4.1 The cosmological density of baryons from 
(D/H) p 

Steigman (2007) has pointed out that in the range of interest 
here the primordial abundance of D is simply related via the 
expression 

10 5 (D/H) P = 2.67(1 ± 0.03) (6/rno) 1 ' 6 (7) 

to 7710 which measures the universal ratio of the densities of 
baryons and photons in units of 10~ 10 : 

7710 s W 10 {n h /n-,) = 273.9 Q. h , h 2 (8) 

The 3% uncertainty in eq. (0 is comparable to the uncer- 
tainties in the nuclear reaction rates used in Big-Bang nu- 
cleosynthesis codes. The conversion from 7710 to Q,b,ah 2 in 
eq. ((§} is accurate to about 0.1% (Steigman 2006). From 
eqs. (6), (7) and (8), we find: 

f7 b ,oh 2 (BBN) = 0.0213 ± 0.0009 ± 0.0004 (9) 

where the error terms reflect the uncertainties in, respec- 
tively, (D/H) p (eq. |6]) and the nuclear reaction rates (eq. [7}. 
Combining the two error terms in quadrature, we have: 

f3 b ,oh 2 (BBN) = 0.0213 ± 0.0010 (10) 

The analysis by Dunkley et al. (2008) of five years of 
observations with the WMAP satellite concluded that, on 
the basis of the WMAP data alone 

«b,ofr 2 (CMB) = 0.02273 ± 0.00062 . (11) 

The two estimates of Qb,oh 2 agree (just) within the errors; 
we also note that the uncertainty from BBN is now com- 
parable to that from the CMB. Given that the best-fitting 
value of fib, oh 2 (CMB) is tied to those of other cosmological 
parameters, it is of interest to consider the effect of includ- 
ing the value of S7b,o^- 2 (BBN) as a prior in the analysis of 
the WMAP 5 data.' 



5 COMBINED CONSTRAINTS ON 

COSMOLOGICAL PARAMETERS FROM 
THE CMB AND (D/H) P 

The five-year WMAP CMB temperature maps and large 
scale polarization maps together provide tight constraints 
on several combinations of cosmological parameters (Hin- 
shaw et al. 2008; Dunkley et al 2008). In order to constrain 
individual parameters, however, it helps to apply external 



data that can break degeneracies. The baryon density, which 
affects the relative heights of the CMB acoustic peaks, is 
partly degenerate with the spectral index of primordial fluc- 
tuations, n s , since WMAP provides precision measurements 
of only two peaks in the temperature power spectrum. In- 
creasing the value of n s increases the height of the second 
peak, which can be compensated by a decrease in the baryon 
density. Hence we can get a better constraint on the spec- 
tral index by combining the CMB data with an independent 
determination of Qb,oh 2 . 

Among the alternative avenues to fib.o which have been 
considered, the primordial abundance of deuterium is the 
most accurate at present. The baryonic acoustic oscilla- 
tions imprinted in the large-scale distribution of galaxies 
currently provide less stringent constraints on Qb,o than ei- 
ther the CMB or (D/H) p and, in any case, really measure 
the combination of parameters fib,o/fim,o/& (e.g. Blake et 
al. 2007). Measurements of other light elements created in 
Big-Bang nucleosynthesis suffer from systematic uncertain- 
ties which are difficult to quantify (in the case of helium) 
or are poorly understood (for lithium), as discussed in re- 
cent reviews by Molaro (2008) and Steigman (2007) (see also 
Simha & Steigman 2008). In considering only deuterium in 
our joint analysis with the CMB we make the implicit as- 
sumption that at present the primordial abundances of 4 He 
and 7 Li are difficult to reconcile with that of D because of as- 
trophysical considerations (such as how to extrapolate from 
measured values in local astrophysical environments to the 
primordial abundances) and do not reflect an underlying 
departure from standard Big-Bang nucleosynthesis. If the 
latter were the case, we would be unjustified in comparing 
the CMB fluctuations with only (D/H) p . 

Parameter constraints from CMB data are usually en- 
coded in a set of samples from the posterior distribution gen- 
erated from the likelihood function by Markov Chain Monte 
Carlo methods. In combining the CMB data with the value 
reached here from (D/H) p , fi b , /i 2 (BBN) = 0.0213±0.0010, 
we assume that the baryon density likelihood is an uncorre- 
cted Gaussian. Since the constraint is only on one parame- 
ter, and consistent with the baryon density inferred from the 
CMB alone, we can use importance sampling to re-weight 
the parameter samples with the extra constraint (see Lewis 
& Bridle 2002 for details). For each sample in the original 
chain supplied by the WMAP team0 we weight the sam- 
ple by exp [-(n hfi h 2 - 0.0213) 2 /(2 x 0.001 2 )]. Marginal- 
ized constraints for individual parameters can be calculated 
easily from the weighted chains; we use the 'GetDist' pro- 
gram provided with CosmoMC (Lewis & Bridle 2002) to do 
this. 

5.1 Standard ACDM cosmological model 

We assume the simplest flat ACDM cosmological model, 
with a power-law purely adiabatic spectrum of linear primor- 
dial curvature perturbations with spectral index n s and am- 
plitude A B at k — 0.002 Mpc -1 , dark matter density Q c fih 2 , 
baryon density Q,b,oh 2 , optical depth for sharp reionization 
r, and cosmological constant density relative to critical Qa,o 

4 Available from: 

http: //lambda, gsfc . nasa.gov / product /map /dr3 / parameters, cfm 



D/H in the most metal-poor DLA 9 




Figure 5. Probability distributions of cosmological parameters deduced from the analysis of WMAP 5 data alone (red dashed line) 
and including f2b,o'i 2 (BBN) as a prior (black continuous line). 
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Figure 6. The points show pairs of values of the baryon density, 
f2b,o^ 2 ! an d the spectral index of primordial fluctuations, n s , im- 
plied by the WMAP data alone, colour-coded according to their 
value of the Hubble parameter Ho (in units of km s — 1 Mpc -1 ). 
The vertical line shows the value of oft 2 (BBN), with the 
shaded regions showing the lcr and 2a bounds from eq. (ID) . The 
contours are the lcr and 2a constraints obtained by combining 
the WMAP data with (D/H) p . 



(with a fiat prior). Marginalized one-dimensional parameter 
constraints obtained in this way are shown in Fig. [5] 

The combined constraint on the baryon density is 
f^b.o^ 2 = 0.0223 ± 0.0005 where the error corresponds to 
the 68% probability from the distribution shown in the top 
left-hand panel of Fig. [S] For comparison, the value obtained 
by combining the WMAP 5 data with the distance mea- 
surements from Type la supernovae and the baryon acous- 
tic oscillations in the large-scale distribution of galaxies is 
^b.o^ 2 = 0.02265 ± 0.0006 (Komatsu et al. 2008). 

Since the (D/H) p measurement prefers values of the 
baryon density towards the lower end of the range allowed by 
the CMB, inclusion of this prior leads to a lower value of the 
spectral index: we obtain n s = 0.956± 0.013, and n s < 0.990 
at 99% confidence (purely statistical errors). Thus, red spec- 
tral tilts are preferred and a Harrison- Zeldovich spectrum 
with n s = 1 is ruled out in the simplest ACDM models. 
Other parameter constraints are almost unchanged, with 
only slight shifts in r, as and Hq towards lower values. 
The combined goodness-of-fit parameter XcS increases by 
about one on including the r2b,o'i- 2 (BBN) constraint, consis- 
tent with expectations for adding one independent param- 
eter. Figure [6] shows the combined constraints on n s and 
&h,oh 2 with and without inclusion of the (D/H) p prior. 

5.2 Constraints on a tensor component to 
primordial fluctuations 

If the spectral index is less than one, some inflationary 
models predict an observable amplitude of primordial grav- 
itational waves. Adding a parameter r (with flat prior) 
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Figure 7. Probability distributions of cosmological parameters deduced from the analysis of: (i) WMAP 5 data alone (red dashed line); 
(ii) WMAP 5 + (D/H) p (blue dot-dash line); (iii) adding to (ii) the constraints imposed by baryon acoustic oscillations in the large-scale 
distribution of galaxies and by the distances to Type la supernovae (black continuous line - see text for relevant references) . 



to measure the ratio of the tensor and scalar power at 
k — 0.002 Mpc -1 , is therefore well motivated (using the re- 
lation between the tensor and scalar tilts for the simplest 
inflation models; see Komatsu et al. 2008). A combination 
of WMAP 5 data and the (D/H) p constraint then gives 
n s < 1.00, and a limit on the tensor amplitude r < 0.26, both 
at 95% confidence. The constraints on these parameters are 
consistent with, and comparable to, those obtained by Ko- 
matsu et al. (2008) by considering WMAP data in conjunc- 
tion with baryon oscillations and supernova distances mea- 
surements, but have completely different systematics. Con- 
sidering all the constraints together - WMAP 5, (D/H) p , 
Type la supernovae (Riess et al. 2004; Astier et al. 2006; 
Wood-Vasey et al. 2007) , and baryon acoustic oscillations in 
the large-scale distribution of galaxies (Percival et al. 2007 ) 
- we arrive at the full joint probability distributions shown 
by the black continuous lines in Fig. [7] We find r < 0.16 at 
95% confidence and n s < 0.994 at 99% confidence. These 
limits strongly constrain possible models of the early uni- 
verse, but are consistent with some of the simplest infla- 
tionary (and other) models. Joint marginalized constraints 
on the various parameters are: O,b,oh 2 = 0.0224 ± 0.0005, 
n cfi h 2 = 0.1130 ± 0.0034, Ov,o = 0.723 ± 0.015, h = 
0.700 ± 0.013, n s = 0.959 ± 0.013, ct 8 = 0.808 ± 0.025, and 
r = 0.818 ± 0.016 {la errors). 



5.3 A closer look at the impact of the errors on 
D/H on the derivation of cosmological 
parameters 

The results presented in Sections 15.11 and 15.21 assumed a 
Gaussian error on f2b,oft. 2 (BBN) derived from the dispersion 
of D/H measurements analysed with the bootstrap method. 
However, values for high-significance confidence limits are 
quite sensitive to the shape of the tails of the distribution 
and it is therefore worth assessing the robustness of our re- 
sults to changes in the statistical model. 

As discussed in Section [3] a Bayesian model that in- 
creases the error bars by some constant and marginalizes 
over the possible values of the constant gives an error es- 
timate on ((D/H) p ) which is comparable to that obtained 
with the bootstrap method. However, the tails of the distri- 
butions have quite different shapes. Fig.[8]shows the result of 
numerically evaluating the distribution of Q^^h 2 , assuming 
Gaussian errors on the log (D/H) measurements, marginaliz- 
ing over the additional error variance with a Jeffrey-like (i.e. 
inverse) prior on its amplitude, and including the 3% error 
(also assumed to be Gaussian) in the relation between rjio 
and (D/H) p (eq.[7]|. The broad tails arise from marginaliza- 
tion over relatively large values of the added error variance, 
and fall off as a power law, rather than exponentially as in 
the bootstrap model. The result is not very sensitive to the 
inclusion or exclusion of the two most outlying data points 
in Fig. H 

Using the marginalized results for additive additional 
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Figure 8. Likelihood distributions for fl b q/i 2 (BBN) for different 
statistical models of the errors in the seven D/H measurements 
considered here (see Table [3$. The solid line is a Gaussian with 
n bfi h 2 = 0.0213 ± 0.0010 (eq. [TOJ. The other curves show the 
result of marginalizing over an increased error variance adding to 
(red dashed) or multiplying (blue dot-dashed) the errors on the in- 
dividual D/H determinations reported in the original works. The 
additional error variance has a Jeffrey-like (regularized inverse 
amplitude) prior. All errors are assumed to be uncorrclated, and 
the results are normalized to peak at 1. 



error variance, the parameter constraints change somewhat 
from the values given in sections 15.11 and 15.21 For exam- 
ple, combining the numerical baryon density likelihood with 
WMAP 5 data for tensor models, the 95% limit on n B 
changes from n s < 1.00 to n s < 1.01 due to the tail of 
higher Qb.oh 2 values that are now allowed. The combined 
tightest constraints change to Q.h,oh 2 = 0.0226 ± 0.0006, 
n s = 0.961 ± 0.014, with n s ^ 1 still just excluded at 99% 
confidence and r < 0.17 at 95% confidence. 

In both the Bayesian and bootstrap model, our results 
depend critically on the assumption that the quoted errors 
on the individual determinations of (D /H) in Table (orig- 
inal and inflated) are statistically independent and Gaus- 
sian. The changes in the values of (or limits On) f^b,0^ ; ^Sj 
and r which we have just noted when changing the statisti- 
cal model of the errors give an indication of the systematic 
error in the result due to the likelihood model, even when 
these assumptions are satisfied. In future, with a larger sam- 
ple of measurements of (D/H) in QSO absorbers, it may be 
possible to test the assumptions that the errors are indepen- 
dent and Gaussian, and perhaps identify the origin of the 
excess scatter in the existing measurements. Such improve- 
ments would allow the primordial abundance of deuterium 
to fulfil its potential in constraining cosmological parame- 
ters. On a different note, the Planck satellite should inde- 
pendently measure the baryon density from the CMB alone 
with ~ 0.6% accuracy in a few years, with measurements of 
several acoustic peaks breaking the degeneracy with n s . 



6 SUMMARY AND CONCLUSIONS 

The lowest metallicity damped Lye* systems are good can- 
didates for the measurement of the primordial abundance of 
deuterium, not only because the gas has suffered little astra- 
tion, but also because they preferentially arise in gas clouds 
with low internal velocity dispersions, facilitating the reso- 
lution of the isotope shift of 82 km s . Furthermore, the 
high column densities involved give detectable D I absorp- 
tion in many lines of the Lyman series. The z^bs ~ 2.61843 
DLA in the spectrum of the bright QSO Q0913+072 is a 
good case in point - we have reported here high S/N ratio 
detections of six D I absorption lines from which we deduce 
log (D/H) = -4.56 ± 0.04. The main contribution to the 
error in log (D /H) is from the uncertainty in the column 
density of Hi, rather than Dl; it may be possible to reduce 
this further in future with more sophisticated modelling of 
the Lyman forest absorption that is superimposed on the 
wide, damped profile of the Lya fine. 

The value of D/H we deduce for this DLA is in good 
agreement with those of four out of the six QSO absorbers 
previously considered to constitute the most reliable set of 
such determinations. The other two differ from the mean of 
the whole sample by ~ 2a, and we have identified possi- 
ble reasons why their errors may have been underestimated 
in the original reports. We propose that the determination 
of the primordial abundance of deuterium is converging to- 
wards the value (log(D/H) p ) = —4.55 ± 0.03, which can 
be considered reliable within the 68% confidence limits. The 
corresponding fi b ,o/i 2 (BBN) = 0.0213± 0.0010 agrees within 
the errors with fi b , ft 2 (CMB) = 0.02273 ± 0.00062. Includ- 
ing the former as a prior in the analysis of the WMAP 5 
data from which the latter is deduced results in a lower 
mean value of the power-law index of primordial fluctu- 
ations, from n s = 0.963 ± 0.015 to n s = 0.956 ± 0.013. 
The effects on other cosmological parameters deduced from 
the analysis of the CMB are more modest. Considering to- 
gether the constraints available from WMAP 5, (D/H) p , 
baryon oscillations in the galaxy distribution, and distances 
to Type la supernovae, we arrive at the current best esti- 
mates tt hfi h 2 = 0.0224 ± 0.0005, n s = 0.959 ± 0.013 (la 
errors), and r < 0.16 (2a limit) for the ratio of the tensor 
and scalar power at k = 0.002 Mpc^ 1 . 

Despite the long integration time, the data presented 
here have not led to a significant change in the estimate 
of (D/H) p compared with the last paper to report new 
observations (O'Meara et al. 2006). This is a good sign; 
together with the agreement between r2b,o^ 2 (BBN) and 
^b,o^ 2 (CMB), it suggests that we have reached a stage 
where the primordial abundance of deuterium and the cos- 
mological density of baryons are sufficiently well determined 
quantities. There is now less urgency to increase further the 
number of D/H measurements in QSO absorbers, although 
a 'prime' sample of only seven data points clearly does not 
allow for complacency. In particular, an expanded sample of 
D/H measures at high redshifts would improve the statisti- 
cal model on which the errors on the derived cosmological 
parameters are based. The most metal-poor DLAs remain 
a fertile ground for studying early episodes of stellar nu- 
cleosynthesis (e.g. Pettini et al. 2008), and it is certainly 
worthwhile continuing to be aware of their potential for fur- 
ther refinements of (D/H) p . 
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APPENDIX A: UVES SPECTRUM OF 
Q0913+072 

For the interested reader, we provide the final reduced and 
co-added UVES spectrum of Q0913+072, as used in the 
work reported here and in Pettini et al. (2008). Relevant 
details are given in Section 2. The spectrum is available 
in its entirety in the electronic edition of the journal. A 
portion is shown here for guidance regarding its form and 
content. 



Wavelength 3, 

(A) 

3.2999990E+03 
3.3000266E+03 
3.3000540E+03 
3.3000815E+03 
3.3001091E+03 
3.3001367E+03 
3.3001641E+03 
3.3001917E+03 
3.3002192E+03 
3.3002468E+03 



Relative Flux 
(Arbitrary Units) 

-4.8096318E+00 
7.2700119E+01 
8.8452042E+01 
6.0643616E+01 
5.9261909E+01 

-8.5954037E+00 
6.8833572E-01 
9.1782742E+00 

-8.3213625E+00 
2.1003305E+01 



la Error 
in Relative Flux 



1.9800955E+01 
2.1541672E+01 
1.9596445E+01 
1.8316561E+01 
2.0739399E+01 
2.0203918E+01 
1.9290850E+01 
1.8983492E+01 
2.0355886E+01 
2.0274757E+01 



a Vacuum heliocentric. 



